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A very simple nonlinear parallel nonautonomous LCR circuit with Chua's diode as its only non- 
linear element, exhibiting a rich variety of dynamical features, is proposed as a variant of the 
simplest nonlinear nonautonomous circuit introduced by Murali, Lakshmanan and Chua(MLC). By 
constructing a two-parameter phase diagram in the {F — lj) plane, corresponding to the forcing am- 
plitude (F) and frequency (ui), we identify, besides the familiar period-doubling scenario to chaos, 
intermittent and quasiperiodic routes to chaos as well as period-adding sequences, Farey sequences, 
and so on. The chaotic dynamics is verified by both experimental as well as computer simulation 
studies including PSPICE. 



I. INTRODUCTION 

Sometime ago, the simplest nonlinear dissipative 
nonautonomous electronic circuit consisting of a forced 
series LCR circuit connected parallely to the Chua's 
diode, which is a nonlinear resistor, was introduced by 
Murali, Lakshmanan and Chua [Murali et al., 1994]. The 
circuit exhibits several interesting dynamical phenomena 
including period doubling bifurcations, chaos and peri- 
odic windows. However, in the paramater regimes in- 
vestigated, it does not exhibit other important features 
such as quasiperiodicity, intermittency, period adding se- 
quences, and so on. It will be quite valuable from non- 
linear dynamics point of view to construct a simple elec- 
tronic circuit which exhibits a wide spectrum of dynam- 
ical phenomena [Lakshmanan & Murali, 1996]. In this 
paper, we wish to point out that a rich variety of phe- 
nomena can be realized with a simple variant of the above 
MLC circut, by connecting the Chua's diode to a forced 
parallel LCR circuit instead of the forced series LCR cir- 
cuit. The resultant system exhibits not only the familar 
period doubling route to chaos and windows but also in- 
termittent and quasiperiodic routes to chaos as well as 
period adding sequences and Farey sequences, to name a 
few. 



II. CIRCUIT REALIZATION OF MLC 
VARIANT CIRCUIT 

The circuit realization of the proposed simple non- 
autonomous circuit, namely, a variant of the standard 
MLC circuit, is shown in Fig. |l|. It consists of a 
capacitor(C), an inductor(i), a resistor(i?), an exter- 




FIG. 1: Circuit realization of the simple non-autonomous 
MLC variant circuit. Here A'^ is the Chua's diode. 
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nal periodic forcing voltage source and only one non- 
linear element(iV), namley, the Chua's diode. In the 
dynamically interesting range, the v — i characteris- 
tic of the Chua's diode is given by the usual three 
segment piecewise-linear function [Chua et al., 1987; 
Kennedy, M.P, 1992; Cruz, J.M & Chua, L.O, 1992]. The 
nonlinear element is added to the familiar forced parallel 
LCR circuit instead of the series LCR of MLC circuit 
[Murali et al., 1994]. The resulting circuit can be consid- 
ered as another important very simple dissipative second 
order nonautonomous nonlinear circuit and a variant of 
the MLC circuit. By applying Kirchhoff's laws to this 
circuit, the governing equations for the voltage v across 
the capacitor C and the current i l through the inductor 
L are represented by the following set of two first-order 
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non-autonomous difFrential equations: 



ih - g{v), 



L 



diL 
dt 



(la) 
(lb) 



where g{.) is a piecewise-linear function defined by Chua 
et al. and is given by g(w) = Gi,v+Q.5{Ga — Gb)[\v + Bp\ — 
\v — BpW which is the mathematical representaion of the 
charactierstic curve of Chua's diode [Kennedy, M.P, 1992; 
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FIG. 2: Two parameter phase diagram in the {F — ui) plane. 

Cruz, J.M & Chua, L.O, 1992]. The values of the various 
circuit parameters chosen for our study are as follows: 
C==10.15 nF, L=445 mH, i?=1745 ohms, Q = l.imKHz, 
Ga = -0.76ms, Gb = -0.41ms, and Bp = l.OV. 

In cquation(l), / is the amplitude and Q is the angular 
frequency of the externl periodic signal. Rescaling Eq. (|^) 
as w = xBp, iL = GyBp, G = 1/R, ui = QG/G, and 
t = tG/G and then redefining r as t, the following set of 
normalized equations are obtained: 



X — F sm{ujt) ^ X — y — g{x), 

- A. 
~ di 



y = f3x, 



(2a) 
(2b) 



where /3 = C/LG^ F = f/Bp, G ^ 1/R. Obviously 
g{x) = bx + 0.5(a - b){\x + 1| - |a; - 1|). Here a = 
Ga/G, b = Gb/G. Note that the two coupled first order 
ordinary differential equation given by Eqs. (^) can also 
be written as a single second order differential equation 
of the Lienard's type in the form 



y + y + Pg{y/I3) + (3y ^ pp sm{ut). 



(3) 



We note at this point chaos via torus breakdown gen- 
erated in a piecewise-linear forced van der Pol equation 
of the form(3) with asymmetric nonlinearity has been 



studied by Inaba and Mori [Inaba, & Mori, 1991] some- 
time ago. However the corresponding circuit uses more 
nonlinear elements than the present circuit. Now the 
dynamics of equation(2) or equivalently (3) depends on 
the parameters a, 6, /?, lo and F. Then for the above 
chosen experimental circuit parameter values, we have 
(3 = 0.05, a = -1.121, b = -0.604 and w = 0.105. We 
use the amplitude F of the external periodic forcing as 
the control parameter and carry out a numerical simu- 
lation of Eqs. (||) either by integrating Eqs. (||) or by 
solving Eq. (^ analytically and numerically, for increas- 
ing values of F. For the above choice of parameters the 
numerical simulation of Eqs. (^ exhibits novel dynamical 
phenomena. 
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FIG. 3: (a) One parameter bifurcation diagram in the {ui — x) 
plane and (b) maximal Lyapunov spectrum (Xmax)- The value 
of the forcing amplitude has been fixed a,t F = 0.28 

To begin with we have carried out an experimental 
study of the dynamics of the circuit given by Fig. |l|. The 
driving amplitude(/) is slowly increased from OV and the 
response of the system is observed to progress through 
a series of transition from periodic motion to aperiodic 
motion. When the driving amplitude / = (correspond- 
ing to the autonomous case), a period- 1 limit cycle is 
observed. By increasing the amplitude from zero up- 
wards, the circuit of Fig. ^ is found to exhibit a quasiperi- 
odic(torus) attractor which then transits to chaos via 
torus breakdown, followed by periodic windows, period- 
adding sequences etc. We have confirmed these results 
also by using the standard 4th-order Rugne-Kutta inte- 
gration routine and carrying out a numerical analysis of 
Eqs. (^ or equivalently Eq. (||) with the rescaled circuit 
parameters of Fig. || as given above. We have observed a 
series of appearence and disappearence of the quasiperi- 
odic (torus) and phase locking(periodic) attractors, and 
period-adding sequences alternatively by varying the am- 
plitude F of the external source at fixed frequency. Illus- 
trated in Fig. ^is the two paramater phase diagram which 
is plotted in the {F — lu) plane. In particular, quasiperi- 
odic motions and period-adding sequences besides the 
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FIG. 4: (a) One parameter bifurcation diagram in the (F — x) 
plane and (b) maximal Lyapunov spectrum (Xmax)- The value 
of the forcing frequency has been fixed at li; = 0.105 



standard bifurcation sequences have been observed in re- 
gions of the lower drive amplitude (F) and frequency {uj) 
values. For example, typical quasiperiodic attractors ex- 
ist in the range F = (0.15, 0.2), and lu =(0.09, 0.16). Sim- 
ilarly, a period adding sequence exists for F = (0.36, 0.38) 
and u! = (0.1,0.115), a period-doubling bifurcation se- 
quence for F = 0.39, and ui = (0.09,0.1), and also type 
I intermittency has been identified for F — 0.38 and 
UJ = (0.10636) . In addition for F = (0.24,0.27) and 
LU = (0.138,0.146), as wch as for F = (0.26,0.29), uj = 
(0.128,0.138) and for F = (0.28,0.33), lj = (0.11,0.125), 
Farey sequences [Kaneko, 1986] exist. The regions of 
chaos are also indicated in Fig. Finally Fig. H(a). 
represents the one-parameter bifurcation in the {uj — x) 
plane, for F = 0.28 which consists of quasiperiodicity, 
chaos, windows, period adding sequences and the familiar 
period doubling bifurcation sequence, intermittency and 
so on. In fig. ^(b) its corresponding maximal Lyaponov 
spectrum in {ui — Xmax) is plotted. Also fig. ^(a) repre- 
sents the one-parameter bifurcation in the (F — x) plane, 
for u; = 0.105, fig. ||(b) depicts the corresponding maxi- 
mal Lyaponov spectrum {F — Xmax) ■ Further, the exper- 
imental results obtained for a choice of circuit parameters 
were also confirmed by using the PSPICE [Roberts, 1997] 
deck available to simulate the behavior of the circuit in 
Fig. 0. In Fig. ^, we have included for comparison the 
chaotic attractor corresponding to F = 0.389^ in Eq. (|l|) 
for experimental and SPICE analysis and corresponding 
numerical simulation for amplitude / = 0.383 in Eq. (0). 





FIG. 5: (a) Chaotic attractor corresponding to F = 0.389V, 
fl = l.llGKHz in Eq. (§): (a) experimental (b) PSPICE and 
(c) corresponding numerical simulation of Eq. (0) 



shown to possess a very rich variety of dynamical phe- 
nomena. In view of the appearance of several ubiquitous 
routes to chaos such as quasiperiodicity, intermittency, 
period doubling, period-adding and Farey sequence in 
such a single but simple circuit, one can make use of the 
circuit in varied investigations on chaotic dynamics and 
applications, including spatio-temporal studies. A more 
detailed analysis of these aspects will be published else- 
where. 



III. CONCLUSION 

In this letter, we have presented a very simple second 
order dissipative nonlinear nonautonomous circuit which 
is a variant of the MLC circuit and carried out. It is 
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